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We reanalyze the prompt muon neutrino flux from gamma-ray bursts (GRBs), at the example 
of the often used reference Waxman-Bahcall GRB flux, in terms of the particle physics involved. 
We first reproduce this reference flux treating synchrotron energy losses of the secondary pions 
explicitly. Then we include additional neutrino production modes, the neutrinos from muon decays, 
the magnetic field effects on all secondary species, and flavor mixing with the current parameter 
uncertainties. We demonstrate that the combination of these effects modifies the shape of the 
original Waxman-Bahcall GRB flux significantly, and changes the normalization by a factor of three 
to four. As a consequence, the gamma-ray burst search strategy of neutrino telescopes may be 
based on the wrong flux shape, and the constraints derived for the GRB neutrino flux, such as the 
baryonic loading, may in fact be already much stronger than anticipated. 



Neutrino telescopes, such as IceCube [T] or 
ANTARES [2], are designed to detect neutrinos 
from astrophysical sources. There are numerous candi- 
date sources, see Ref. [3] for a review and Ref. [I] for 
the general theory. We focus on the prompt emission 
of gamma-ray bursts (GRBs) in this letter, where 
photohadronic interactions are expected to lead to a 
significant flux of neutrinos [5]. So far, no extraterrestrial 
high energy neutrino flux has been detected yet. That 
is, for sources optically thin to neutrons, consistent with 
generic bounds O [7] which are just being touched by 
tceCube. The search for GRB neutrinos has been driven 
by analytical estimates for the shape and normalization, 
the simplest one being the Waxman-Bahcall (WB) 
flux [B]. More recent analyses, such as the stacking 
analysis in Ref. [8], relating the neutrino flux to the 
observed gamma-ray flux, are based on the analytical 
generalization of this flux for arbitrary input param- 
eters following Ref. [S]. These calculations typically 
approximate the A(1232) resonance for the charged pion 
production 
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in some form. However, the GRB neutrino flux compu- 
tation has been updated over the last ten years from the 
particle physics point of view by improving the descrip- 
tion of the photo-meson production processes, and it has 
been obvious there is a substantial impact from mag- 
netic field effects and flavor mixing on the neutrino flux 
as well; see, e.g., Refs. [THl - TH] . In this letter, we make 
the impact of these effects very explicit by revising the 
often used WB reference flux from Ref. [6]. We include 
the relevant pion production modes and neutrinos from 
kaon and neutron decays. We treat the magnetic field 
effects on each charged particle species explicitly, and we 
include flavor effects/flavor mixing. Note that we keep 
our considerations as independent of the astrophysical 
source model as possible to factor out the particle physics 
effects, which are much better known than the details of 
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FIG. 1: The WB flux from Eq. Q (thin dashed curve), 
the numerically reproduced flux using the A + resonance only 
(lower solid curve), and the WB flux including higher reso- 
nances, direct production/t-channel processes, and multi pion 
production (high energy processes), which are successively 
switched on, leading to the final upper solid curve. Here the 
Vfj, flux from n + and ir~ decays is considered. The normal- 
ization of our result to the numerically reproduced WB flux 
(gray dashed curve) is described in the main text. 



the astrophysical model. The purpose of this letter is to 
demonstrate how the original WB flux changes in both 
shape and normalization effect by effect, and where the 
main impact comes from. We also discuss the impact on 
data analyses. The technology used in this letter is based 
on Refs. fl4l 1151. where details can be found. 



In the standard picture, protons collide with photons, 
possibly from synchrotron emission of co-accelerated elec- 
trons or positrons (see, e.g., Ref. [H]), leading to pion 
production by processes as, for instance, Eq. Q. The 
charged pions then decay further into neutrinos, such as 
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of the WB flux, consider only the from pion decays 
for the moment. It is often assumed that the target pho- 
ton field corresponds to the observed prompt GRB flux, 
which is typically parameterized by dN 1 (E) / dE oc E a i 

for E < £ 7 ,brcak and dN 1 (E) / dE OC E@~t for E > £ 7 ,break 

in the observer's frame, where a 7 ~ —1, /3 7 ~ —2, and 
the break e 7i broak at a few hundred keV. If the protons 
are injected with a power law with injection index two, 
one obtains for the prompt GRB neutrino flux, referred 
to as "WB flux" , 
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with a v = —(If — 1 
£?, ~ 10 5 GeV and si 
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+1, P u = —®i — 1 — 0, 
10 7 GeV. For the analytical 



estimates of the break energies, we follow the treatment 
in Ref. [5], assuming that T = 10 2 5 and z = 2; see, e.g., 
Refs. E]- The first break energy e° v can be related 
to 

£7, break from the threshold of the photohadronic inter- 
actions at the source. As a minor difference to Ref. [9], 
where heads-on collisions between photons and protons 
are assumed for the threshold, we include the effect that 
the pion production efficiency peaks at higher center-of- 
mass energies (see Fig. 4 in Ref. [TOJ) to match our numer- 
ical results. This leads to a factor of two higher photon 
energy break (14.8 keV) in the source frame to match the 
e v ~ 10 5 GeV for the chosen parameter set. The second 
break comes from pion cooling in the magnetic field. It 
can be computed from the energy where the pion decay 
rate equals the synchrotron loss rate. In order to repro- 
duce e s v ~ 10 7 GeV, one has B ~ 3 • 10 5 G. Note that, in 
the light of recent Fermi data, it is not clear how "typi- 
cal" this parameter set is, which, however, does not affect 
the logic of this letter. As another relevant parameter, 
we choose the maximum proton energy by balancing syn- 
chrotron loss and acceleration rates with an acceleration 
efficiency of 10% [TO]. For the expected normalization of 
the flux in Eq. Q, we use [5] (updated in Ref. [TO] ) 

■GeVcm^s^sr- 1 (3) 
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per neutrino species (y e , v^, or D^). After flavor mix- 
ing, the combined muon neutrino and antineutrino flux 
is, again, approximately given by Eq. ^ [2U]. This esti- 
mate is based on the assumption that GRBs are a domi- 
nant cosmic ray source in which the high energy protons 
dissipate a fraction / w < 1 of energy into pion production 
before leaving the source. If no neutrino flux at the level 
of Eq. ([3]) is observed, it means that effectively the prod- 
uct of f n and the fraction of energy in protons (baryonic 
loading) becomes stronger constrained [9] - and there- 
fore the hypothesis of GRBs being the dominant cosmic 
ray source. We choose /„. = 0.2 for the following figures. 
For the numerical treatment of the photohadronic inter- 
actions, we follow Ref. [TO] (Sim-B), based on the physics 



of SOPHIA [TO] and the weak decays in Ref. [TO], includ- 
ing the helicity dependence of the muon decays. The 
energy losses and other production modes are treated as 
described in Ref. [TO] . We assume that synchrotron losses 
are the leading energy loss mechanism, which means that 
only the product of the proton and photon densities is re- 
quired, see Eq. (7) of Ref. [13] . Therefore, our results are 
independent of the baryonic loading and GRB model de- 
tails. We also assume that the source is optically thin 
to neutrons, which means that secondary interactions 
are neglected. There are limitations to these assump- 
tions, such as if the protons cool significantly by photo- 
hadronic interactions (see, e.g., discussion in Ref. |11|). 
However, such processes cannot be included in a model- 
independent way, because they require separate knowl- 
edge on the proton and photon densities. We show in 
Fig. [l] the WB reference flux from Eq. ([2| as thin dashed 
curve for from charged pion decays only. In order to 
normalize our flux to this curve, we need to take into 
account that the assumptions for the A resonance vary 
in the literature. Therefore, we first of all reproduce 
the WB flux numerically by including the synchrotron 
cooling of the pions explicitely, leading to the thick gray 
dashed curve "WB A + -approx." . Here we use the same 
cross sections, pion multiplicities, and inelasticities as in 
Ref. [5], and we choose the normalization such that the 
energy going into neutrinos is the same as for the analyt- 
ical estimate. Note that for this curve, the second break 
is automatically reproduced by magnetic field effects and 
not put in by hand, which results in a small pile-up effect 
at the plateau. By this choice, the product of proton and 
photon density normalizations is fixed, and we can use 
the input spectra to compute the effects of the more re- 
fined interaction model. We show this by the solid curves, 
where higher resonances, direct (t-channel) production, 
and multi pion production are successively added to the 
actual A + (1232) resonance process in Eq. Q; see, e.g., 
Refs. [TO] [TO]. The final result exceeds the WB estimate 
by a factor of a few, especially at high energies. The ad- 
ditional tilt of the spectrum comes from the multi pion 
cross section staying approximately constant for high in- 
teraction energies. In addition, note that all processes 
other than the actual A resonance include 7r~ production 
and two or multi pion production modes. From Fig. [l] 
one can read off that the WB approximation basically 
includes the effect from direct production at low ener- 
gies and higher resonances at high energies, but the high 
energy (or two and multi pion) contributions are clearly 
underestimated. In addition, one can read off that the 
effect of the additional production processes can lead to 
an up to one order of magnitude change of the flux, de- 
pending on the assumptions on the A resonance in the 
literature. 

Apart from pion decays, neutrinos are produced from 
muon decays in the pion decay chain. In addition, kaons 
produced in the photohadronic interactions similar to pi- 
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FIG. 2: Total electron neutrino (left panel) and muon neutrino (right panel) flux before flavor mixing (thick solid curves), 
where the neutrino and antineutrino fluxes are added. The individual contributions to these neutrino fluxes from tt, fi, n, and 
K decays are shown as well. 



ons may decay into neutrinos |21j . The main qualitative 
difference among charged pions, muons, and kaons are 
their different masses and lifetimes, leading to different 
energies of the second (synchrotron) break in Eq. ([2]); see, 
e.g., Fig. 3 in Ref. [15]. This effect has interesting impli- 
cations for the flavor ratio of the neutrinos, which changes 
as a function of energy |T2] ; see also Ref. [55] . Finally, any 
neutrino flux will be accompanied by a neutron flux, as it 
is obvious from Eq. ([I]). These neutrons are, however, not 
stable. If the neutrons do not interact, they will decay 
either within or outside the source by n — > p + e~ + D e , 
leading to cosmic rays and (inevitably) to an additional 
(almost coincident) v e neutrino flux. We show the total 
electron neutrino (left panel) and muon neutrino (right 
panel) flux before flavor mixing in Fig. [2j where the neu- 
trino and antineutrino fluxes are added. The individual 
contributions to these neutrino fluxes from tt, fi, n, and 
K + decays (we only consider the leading kaon contribu- 
tion mode) are shown as well. The WB flux from Eq. ([3| 
is given as reference for the corresponding number of neu- 
trino species. In the right panel (muon neutrinos), one 
can clearly see the hierarchy in the second break energy 
among neutrinos from /i, tt, and K decays. In the left 
panel (electron neutrinos), the main contribution comes 
from muon decays. However, neutron decays show up at 
low energies. 

In order to obtain the final muon neutrino flux relevant 
for muon tracks in neutrino telescopes, the total elec- 
tron and muon neutrino fluxes in Fig.[2]are superimposed 
by flavor mixing (averaged neutrino oscillations) [20], see 
Fig. [3] The shaded band indicates the 3a allowed range 
of the total flux from current mixing parameter uncer- 
tainties [23]. In fact, if the combined knowledge from 
the Double Chooz, Daya Bay, T2K, and NO^A is ap- 



plied [53], as expected in about 2015, this band becomes 
hardly visible anymore. Therefore, mixing parameter un- 
certainties are less relevant for the GRB analysis, espe- 
cially at the lower break (unless flavor ratios are con- 
sidered). Comparing the final result (upper thick solid 
curve) with the WB flux (thin dashed curve), we notice 
that the expected neutrino flux is about a factor of three 
to four larger than the WB flux at 1 PeV. In addition, 
we show a rescaled version of the final result (thin solid 
curve) to illustrate the impact on the spectral shape com- 
pared to the WB flux. It is obvious from this compari- 
son that the shape of Eq. ^ cannot be used for realistic 
data analyses or to search for point source GRBs. For 
instance, the first break, to which AMANDA and Ice- 
Cube are most sensitive to, has basically disappeared in 
its original form. On the other hand, magnetic field and 
flavor effects lead to a characteristic double peak struc- 
ture, one could search for if a few bursts dominated. In 
addition, note the high energy excess coming from kaon 
decays, which increases the flux by at least one order of 
magnitude. At about 10 8-9 GeV, horizontal air shower 
experiments, such as Auger [25], in fact have the best 
sensitivity, where the flux shown in Fig. [3j is representa- 
tive for or v T events. Although the full-scale diffuse 
flux sensitivity is considerably above the expected neu- 
trino flux, Auger may detect a flux from GRB kaon de- 
cays especially if a bright burst out-shines the cosmogenic 
neutrino flux. Therefore, neutrino point source studies 
should be initiated by these experiments. 

In summary, we have revised the WB neutrino flux, of- 
ten used as a reference GRB flux, by including the most 
relevant neutrino production processes, and by treat- 
ing magnetic field and flavor effects explicitly. We have 
used as few assumptions as possible on the astrophysi- 
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FIG. 3: Total muon neutrino flux after flavor mixing (dark 
thick solid curve). The individual contributions to this flux 
from muon and electron neutrinos are shown as well (from 
thick curves in Fig. [2). The WB flux from Eq. Q is shown for 
reference, corrected by flavor mixing. In addition, a rescaled 
total flux is shown to illustrate the impact on the spectral 
shape. The shaded band shows the 3<r allowed range of the 
total flux from current mixing parameter uncertainties [23] . 
Here also the 10 year (extrapolated) full-scale limits from 
IceCube [T] (for an E~ 2 diffuse flux above the atmospheric 
neutrino background) and Auger [25] (differential limit) are 
estimated at the 90% CL. 



cal source model. We have demonstrated that the flux 
normalization increases by a factor of three to four with 
respect to the initial assumptions, and that the spec- 
tral shape exhibits a double peak structure qualitatively 
different from the WB flux. The main impact are ad- 
ditional neutrino production modes and magnetic field 
effects, which act differently on the charged secondary 
particle species. The revised spectral shape may allow 
for new search strategies for GRB neutrino fluxes. Since 
current state-of-the-art multi-messenger stacking analy- 
ses, such as Ref. [5], where the expected neutrino flux is 
computed from the observed gamma-ray flux on an event- 
by-event basis, rely on the assumptions on the photo- 
hadronic interactions, the non-observation of a flux may 
have stronger constraints on the energy equipartition be- 
tween protons and electrons than anticipated. 
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